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ABSTRACT 

Light refracted by the planet's atmosphere is usually ignored in analysis of 
planetary transits. Here we show that refraction can add shoulders to the transit 
light curve, i.e., an increase in the observed flux, mostly just before and after 
transit. During transit, light may be refracted away from the observer. Therefore, 
even completely transparent planets will display a very similar signal to that 
of a standard transit, i.e., that of an opaque planet. We provide analytical 
expression for the amount of additional light deflected towards the observer before 
the transit, and show that the effect may be as large as 10 -4 of the stellar 
light and therefore measurable by current instruments. By observing this effect 
we can directly measure the scale height of the planet's atmosphere. We also 
consider the attenuation of starlight in the planetary atmosphere due to Rayleigh 
scattering and discuss the conditions under which the atmospheric lensing effect 
is most prominent. We show that, for planets on orbital periods larger than 
about 70 days, the size of the transit is determined by refraction effects, and not 
by absorption within the planet. 

Subject headings: planetary systems - techniques: photometric - planets and 
satellites: general 



Introduction 



While the maj o rity of extrasolar p lanets were found using the radial velocity technique 
(IMayor et al.l 120041 ; iMarcy et al.l 120081). the method of finding planets by observing their 
transit in front of a star (jBorucki fc Summers! 1 19841 ) is now becoming a powerful tool. The 
transit method is based on the observation of a small drop in the brightness of a star, that 
occurs when the orbit of a planet passes ('transits') in front of it. The amount of light lost 
depends on the ratio of the radius of the planet to the radius of the star. It is of order 0.01% 
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(for earth-size planets) and 1% (for Jupiter-size planets). Transiting extrasolar planets reveal 
more information. Since the star's size can be estimated from spectroscopic observations, 
the planet's size ca n be determi ned. The main effect of planetary transit can be modeled 
geometrically (e.g., ISeagerl 120081 ). Both star and planets are projected on the plane of the 



sky, producing two circles of radius R* and R p . Then, the fractional overlap between the 
two determines the fractional decrease in stellar light. Within transit, the size of the effect 



is (R„/R*) 2 . The high a ccuracy obtained, especially by space instruments (IBrown et al. 



2001 



Borucki et al.ll2009l ) demand more advanced models, which take into account the non 



uniform stellar brightness over the stellar disk. This results in a curved, rather than a flat, 
lightcurve around the eclipse minimum. In addition, reflection and thermal emission from 
the planet (|L6pez- Morales & Seagerll2007l) ca us e an out of eclipse sin usoidal variation in the 



flux (IBorucki et al. 



2009; 



Knutson et al. 



20091 ) . iBurrows et al.l (120071 ) stress that care should 



be taken in the definition of the planet radius. While the canonical definition uses radial 
optical depth of r = 2/3, the light in a transit is absorbed tangentially, and therefore passes 
a longer path through the atmosphere of the planet than a radial ray would if it were to 
penetrate to the same depth. This results in an increase of the apparent radius of the planet 
by about 5 planetary scale heights. 



Seager fc Sasselovl (120001 ) showed that the lengthening of the light path due to refrac- 



tion in the atmosph ere is small and therefore does not change the absorption significantly. 
Hui fc Seager! (120021 ) considered refraction by planetary atmospheres, taking the effects of 
oblateness and absorption into account. In this paper we consider the effects of refraction 
in the planetary atmosphere independent of absorption. Outside of the transit, refraction 
may redirect light through the planetary atmosphere towards the observer. This results in 
increased observed brightness especially right before and after a transit. In contrast, during 
transit, refraction deflects light that could otherwise arrive to the observer. This results in a 
diminishing of light even for a completely transparent planet, where no light is absorbed. We 
show that for planets occupying obits with large enough semimajor axis, this refraction effect 
during transit results in an effective planetary radius, Rp, that is larger than that dictated 
by opacity. Therefore, in some cases we would detect transits of completely transparent 
planets. 

The structure of the paper is as follows. In §2] we lay out our model and simplifying as- 
sumptions. In £j3]we derive the deflection angle through a spherically symmetric exponential 
atmosphere along with some analytic results for light trajectory in such systems. Section §4] 
discusses the enhancement of light outside transits and the decrease in light during transit. 
In §S] we take the effects of scattering into account. We calculate the condition where scat- 
tering eliminates the additional light before and after transits and where it dominates the 
effective size of the planet during transits. 
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2. Model assumptions 

We assume a planet in hydrostatic equilibrium with a spherically symmetric exponen- 
tially decaying density profile. Since the deviation of the index of refraction of gas from unity 
5n = n — 1 is proportional to its density, we also have an exponentially decaying 5n. The 
index of refraction n as function of the distance from the center of the planet, r, is therefore: 

/ r - R x \ 

n{r) = 1 + 5n(r) = 1 + exp I — — j , (1) 

where H is the scale height of the planetary atmosphere. Here we used the parameter Ri as a 
nominal normalization of the index of refraction. It is the distance from the center where 5n 
is unity. In planets like Jupiter, it is roughly halfway to the center of the planet. Therefore, 
Ri is many scale heights inward of the usual definition of a planet's radius (e.g., where the 
optical depth radially outward is r = 2/3). The relation between R\ and the radius of the 
planet will be clarified in §4.21 and in the summary. The scale height H is related to the 
planet's atmospheric temperature T p and its surface gravity g by H = where fee is the 
Boltzman constant, T p is the planet's atmospheric temperature, and /i is its mean molecular 
weight. Note that we have used an exponential profile throughout the planet, but only the 
outer atmosphere matters. 

We are mainly interested in light trajectories that pass a few scale heights near the 
radius of the planet, meaning in its atmosphere. Inside the atmosphere, the light is diverted 
due to the changing index of refraction. Our analytical results will be facilitated in the limit 
that the scale height of the planet is much smaller than the planet's radius which in turn is 
much smaller than the stellar radius i?*: H <C R p <C R*. Finally, we assume that the latter 
is much smaller than the semimajor axis of the planet a. This last assumption implies that 
the angular deflection of the light, that bends over the planet's atmosphere and reaches a 
viewer at infinity isA6~i?*/a^Cl. 



3. Light trajectories in a spherically symmetric medium 

3.1. Order of magnitude 

For small deflections, the innermost scale height towards which the light was directed 
dominates the deflection. The distance traversed by an unbent ray within that scale height 
is 2y 'R p H. By definition of a scale height, or equation ([1]), 5n changes significantly over 
a radial distance H. Therefore, in a wide light ray, the part traveling at the top of the 
scale height will be ahead of that at the bottom of the scale height by a distance of order 
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(n — 1) «/ HRp. We can therefore estimate the deflection angle to be (n — 1)^J R p / H . This 
calculation is valid for deflections which would shift the ray less than a single scale height 
over a distance of ^R P H. Therefore, for this result to hold, the deflection angle must be 
less than 50 <C y/H/ ' R p . Such small deflections can direct light towards the observer before 
transit only if R*/a <C a/ H/ R p . 



3.2. Fermat's Principle 

An elegant way of calculating the trajectory of light in the medium is by using Fermat's 
principle. Following the spherical symmetry of the refraction index the trajectory takes place 
on a plane. Hence, we may use polar coordinates: 

n(r) 



dt = -^Vrfr 2 + rW (2) 

c 

and Fermat's principle reads 



*M = */^Jl + r»(§)V = 0. (3) 



Using the Euler-Lagrange equation, we obtain 

d n(r)r 2 f 



0. (4) 



or 

- = T b (5) 

dr y/r*n 2 {r) - b 2 r 2 

Here b is a constant of integration set by the initial conditions. Far away from the planet, 
n ( r )r^oo — 1; and the light travels on a straight line. Solving for n(r) = 1 gives a straight 
line trajectory with a minimal radius b and therefore b is the impact parameter. Taking 
the limit of equation (JHJ) at large r we obtain d6/dr = b/r 2 , and therefore b is the impact 
parameter of the straight line that describes the trajectory far away from the planet, i.e. the 
impact parameter of the non deflected incoming light ray. 

For a non-homogenous index of refraction, (n ^ 1), the impact parameter b and the 
minimum radius r min do not coincide. Since d6/dr diverges at the minimum radius, we can 
find the minimum distance by setting the denominator in equation to zero. We obtain 

^mm^(^mm) = b (6) 
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An alternative way to describe the light trajectory, is to denote by a (0 < a < 2tt) 
the angle between the propagation direction of the light ray and the radial direction. Then, 
tana = rd9/dr. Therefore, equation (jSj) can be elegantly written as 



Equation (jfJJ) is then a special case of equation (jTJ) that describes the turning points, where 
the light ray is tangential so, by definition, a = tt/2. The functional form of rn(r) therefore 
determines the properties of the light trajectories. Interestingly, for our assumed density 
profile, the function rn(r) increases linearly steeply, reaching a maximum at about f = H, 
then decays roughly exponentially, reaching a minimum around f = R± + H In Rx/H, where 
fn(r) = r + H . Outward of that, it increases roughly linearly with rn(r) = r. Note that 
we have assumed that n — 1 grows exponentially all the way to the center. Realistically, 
the density profile around the center grows more slowly. Then, the maximum of rn is not 
at f = H, but much farther out. The minimum we found at f = Ri + H In Rx/H is up 
in the atmosphere, where the exponential assumption is realistic, and therefore is correctly 
estimated. 

We therefore conclude, that a light coming from infinity, with impact parameter b, will 
be deflected from the atmosphere of the planet, if b £ f+H. In contrast, for b ^ f + H light 
rays will deeply penetrate the planet and, may only be reflected after arriving at r r = H. 
We note that all such rays will be absorbed in the planet, under any realistic opacity. We 
now briefly mention two other solutions to equation (JSJ) which are fully contained within the 
planet, and do not emerge out of it. 



a light ray can have a circular orbit around the planet. However, such orbits are unstable, 
and will ultimately either escape the atmosphere or penetrate to large depth. 



rn(r) = b sin a. 



(7) 



3.2.1. Circular Trajectories 



At the minimum of rn or at 



r = r"S R x + H In Rx/H 




3.2.2. Trapped Rays 



Since rn may never increase above b (as sin a < 1) we can have trapped rays that 
propagate between two radii, one close to the atmosphere, and one close to the center of the 
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planet. Such rays will be trapped between r 1 and r 2 satisfying: 

T < r\ < r < r2 < r (9) 
Again, these are likely to be absorbed for any realistic planet. 

3.3. The limit of small deflections 

We determine the angular deflection of light traveling through our spherically symmetric 
medium by integrating equation (jSJ). For convenience, we subtract the result of a light 
traveling within a homogenous medium, n — 1 with the same impact parameter b: 

A0 = 2 /"(*£=- (10) 
Jr min ^Jr 2 n 2 (r)-b 2 ^/r 2 -b 2 

Now, using a change of variables, x = and using flSI) we get 

f min 

A0 = 2 r°( . 1/X = - } IX )dx (11) 
h Jx 2 ^^ - 1 Vx 2 -1 

V ™(»'min) 2 

Under the assumption that r min / H S> 1 yet Snr min /H C 1 we can expand ^ to the 
first order in Sn and get an expression for the angular shift: 



AQ = 6n(b)^2n^ (12) 

This is the basic result we will use next. It is equivalent to what we found using order of 
magnitude derivation in §3.1[ but includes an additional factor of \/2~7r. 



4. Lensed Stellar Images 

Equipped with an analytic expression for the angular deflection, we determine below 
the star's images to a viewer at infinity. Using these images we then estimate the change in 
the observed luminosity due to atmospheric lensing. 

A light arriving from any point on the star could pass by both sides of the planet and be 
deflected towards the observer at infinity. For example in figure [H point A has two images 
A' and A". If the center of the planet is in between the source and the image as projected on 
the sky we refer to the deflection as deflection through the far side of the planet. Otherwise, 
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if, in projection, the source and the image are both on the same side of the center of the 
planet we refer to the deflection as deflection through the near side of the planet. Out of 
transit, the deflection through the far side of the planet would yield a secondary image of the 
star that appears within the planet's atmosphere and has a crescent shape. The light that is 
deflected though the near side of the planet and reaches an observer at infinity, is generally 
way up in planet's atmosphere and is hardly deflected given the exponential nature of the 
index of refraction (for example, points A" and C" in figure [I]). These rays produce the 
simple, almost unperturbed, disk-like image of the star away from transit. The secondary, 
crescent like, image changes the ordinary light curve away from, but close to, transits. In 
order to estimate the change in the observed flux we find the area of the crescent image as 
refraction preserves the specific intensity. 



We denote the star radius as R+ } the semi-major axis of the planetary orbit as a and 
the projected distance between the star and the planet on the major-axis as XR*. 

Now we must formulate a proper way of integrating over the star's edge which will be 
mapped, due to refraction through the far side of the planet's atmosphere, to the crescent 
shaped image. Using Green's theorem for a closed curve integral we can calculate the crescent 
area. If we denote the crescent's area and the distance of each point on the crescent from 
the planet's center by S' and b respectively, then 



4.1. Evaluation of secondary image size 




(13) 



To evaluate the integral, we change integration variable from a' to a using 



tan a' = 



sin a 



(14) 



X — cos a 



For the light to arrive at the observer, the deflection angle must satisfy 



sin A0 



l + b 



(15) 



a 



where from Figure we see that 




(16) 



We can therefore obtain b(a) by equating equations f|T2|) and f fl5|) . While in §4.31 we nu- 
merically solve these equations for b, analytical expressions are given in §4.41 and §4.51 under 
some limits. 
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Fig. 1. — A schematic diagram showing the crescent shape of the secondary stellar image. 
The upper panel shows the projection on the sky while the lower panel shows the plane that 
contains the observer and the center of the planet and star, allowing to follow ray trajectories. 
The points A, B, C and D are mapped by the atmospheric lensing onto points A' B' C and 
D' . On the lower panel, points A and C are also shown to arrive practically undeflected to 
points A" and C" . An arbitrary source point in the star located at an angle a, is mapped 
into a' in the atmosphere of the planet. 

4.2. The Effective Radius of The Planet 

We have used the parameter Ri to normalize our density or index of refraction profiles. 
During the transit itself, light is deflected away from the observer. We can now calculate the 
effective radius of the planet as the radius of the dark spot in the star's image caused by this 
deflection. In the middle of a central transit, i.e., when the center of the planet covers the 
center of the star as projected on the sky, the deflection angle that determines the effective 
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radius is given by ABct = (R* + Rp)/a- This is the deflection of light coming from the edge 
of the star and refracted in the far side of the planet in order to reach the distant observer. 
Using this deflection angle we obtain 

R p ee b (Ae CT ) - R x + | In (^r) ■ (17) 

This is an implicit equation since R p appears on the right hand side as well. However, since 
it only appears in the logarithm, any value for the radius can be used there with negligible 
errors. Exact solutions can be found iteratively with ease. Note, that the above is not a 
standard definition of the planet radius as it depends explicitly on its distance from the star. 
Yet, our R p as defined above is quite similar to the planet's radius defined in other ways. 
Comparison between our effective radius for occultations R p with other definitions of the 
radius of the planet related to atmospheric absorption or scattering are given in §5j This 
definition of the planet radius, allows for a simpler equation of the deflection angle at small 
deflections: 

Ae ^exp(^) (18, 

4.3. Images by numerical calculation 

We numerically integrate equation ( !T3|) and plot the images of the star as seen by an 
observer at infinity. In order that the effect will be easily noticed we chose the following (not 
realistic) parameters a = 100i?*, R* = 8Ri and R\ = AOH. For these parameters, equation 
f fT7|) then implies R p = 1.7 Ri. The image shapes are plotted in Figures (jSj) and (j3J). We can 
see that as X diminishes the secondary image becomes larger. 




Fig. 2. — The image of the star as refracted through an exponential planetary atmosphere 
plotted for different values of planet star separation, X. The axis are measures in units of 
star radius. For X > 1, an additional crescent shape image is seen through the planet's 
atmosphere. When X < 1 the star is being eclipsed even if the planet is completely trans- 
parent. 
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Fig. 3. — The crescent as it appears in the planetary atmosphere for different values of 
the separation parameter X. In this figure the origin is fixed to the planet center, and the 
coordinates are in units of the planet's radius R p . As X decreases the crescent moves away 
from the planet, becomes wider and covers a bigger angle. 

For X < 1 the planet occults the star and an eclipse will be observed. To understand 
that, we take as an example the case where the center of the planet and the center of the 
star coincide in projection on the sky. In this case, every point on the limb of the star will 
either be almost unperturbed when it is deflected from the near side of the planet, while 
deflections from the far side will pass at a distance of R p from the center of the planet (per 
the definition of R p ). Therefore, a dark spot will be observed in the center of the star, with 
radius R p . The change in observed flux affected by the eclipse is given by the ratio of the 
areas S'/S = R p /R 2 . Symmetry, of course, dictates a bright spot at the center of the dark 
one since a light ray that would come from the center of the star passing through the center 
of the planet will arrive without deflection to the observer. However, as we have discussed 
below equation © such rays penetrate so deep into the planet and are always absorbed. We 
therefore ignore them. 

This decrease of (R p / R*) 2 for transparent planets is similar to the one of opaque planets. 
Transparent planets produce similar transits to opaque ones! The light curve that correspond 
to the parameters mentioned earlier is shown in Figure (jl]). The crescent shaped secondary 
image of the star does not yield an analytical expression for a generic of X. An analytical 
expression for the crescent's size can be obtained under several limits. First we study the 
case where the planet is far away from occultation regime, which corresponds to X ^> 1, and 
then we study the crescent size just before transit, where R p / R* < X — 1 -C 1. 
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Planetary transit light curve 




X in units of star radius 



Fig. 4. — The light curve of a star during a completely transparent planetary transit while 
considering atmospheric lensing, for a = 100i?*, R* = 8R P and R p = 100H. Just before 
and just after the transit, the stellar light is increased due to refraction in the atmosphere 
of the planet. Analytical expressions for these shoulders are given in equations ( )23l) . ( 1241) & 
( 130]) . The curved bottom of the transit appe ars without tak i ng lim b darkening into account. 
Similar curvature in transit was obtained by lHui fc Seagerl (120021 ). 



4.4. Away from occultation regime 

From Figure (CQ), we obtain 

. R± r—z — ■— XR± / cos a 

A6 = — VX 2 - 2Xcosa + 1 = 1 — 

a a \ X 



(19) 



We also have an expression for A0 from ffT8]) . Since 6 ^> H and both Rp and b differ by a 
few .ff's we solve for 6(a) as 



R p - H In (X - cos a) = R p - H In X + — cos a 



X 



Noting that for (X ^> 1) equation (fT4"l) can be approximated as 



. cos a 
da = — ——da 



We can estimate the secondary image size 

S' = 



X 



b 2 , . TiHR p 
2 da= ^ 



(20) 



(21) 



(22) 



Now if we assume that the light is not scattered or absorbed, and ignore limb darkening of 
the star, the brightness will be uniform and identical over both the direct and lensed stellar 
images and the fractional increase in light arriving to the observer due to atmospheric lensing 
is 



S' H R p 
J = X 2 R 2 J 



(23) 
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or, substituting the expression for H, 



,S" / n \ -1/2 / R x 3 

— = 6.4 x 10 



5 VlAU/ \Rj 



4.5. Near occultation regime 

Observing the numerical plot of the planetary transit light curve shown at Figure (j4]), 
and the crescent images plotted in Figure (j2j) and ([3]) , one can easily notice that the effect is 
most substantial when X is near unity. Hence, we derive now asymptotic expressions in the 
limit of j^- < X — 1< 1. In this limit the crescent radial width does not change significantly 
until its edges at a' of almost ±~. The area of the crescent is therefore irHRpAb where Ab 
is it's central width. To estimate Ab, we use equation (fTBj) : 

b^ Rp- H In (Ae-^-j (25) 



Following the notation of Figure (Op), we get 

= 771n I 

AQ C > 



Ab^b A ,-b c , = Hln[^^\ (26) 



Noting that 

AQ A , = v ^ ' ' ^ p A6c = — ' ^ P (27) 



(X + 1) R* + Rp A ^ (X - 1) + R p 



together with < X — 1 <C 1 and R+ ^> R p we finally get 



The crescenst's area is then 



' T/V/]n( ^_ i _^ ) (29) 



so that the fractional increase in observed flux is given by 

S' R n H I 2 \ . 

^r = ^ln jr- . 30 

S Rl \X-1 + ^J 

This is up to a factor of In « 2.3 larger than what the extrapolation of equation (20) to 
X = 1 would imply. Figure (jSJ) shows how the asymptotic estimates fl23|) and fl30|) coincide 
with the numerical calculation. 
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tL \ 2 3 4 5 

X-separation in units of stellar radius 



Fig. 5. — The shoulders in transit light curve, ignoring extinction. The numerical curve 
is plotted as blue line for a = 100R+, R± = 10R P and R p = 100H. For comparison, our 
analytical estimates far away from transits (green line, equation (1231) ) and very close to 
transit (red line, equation ( 13"U|) ) are also shown for the same parameters. It can be seen that 
these provide excellent approximations within their range of validity. 

4.6. Conservation of light? 

We have shown so far that light is added to the observer when the planet is out of 
transit. There, the effect is of order HR P /R 2 , and when averaged^ over all possible positions 
of the planet around the star we have HR p /a 2 . However, during transit, refraction causes 
a decrease of order (R p /R*) 2 . Averaged over all possible planet positions we have R 2 /a 2 . 
The two effects, therefore, do not cancel. To settle this apparent discrepancy, one must take 
into account the light that penetrated deeply into the planet. The amount of that light is 
R p /a 2 , but it is being spread roughly equally over 4n steradians as viewed from the planet. 
Therefore, for a non realistic, truly transparent planet, the planet will shine out of transit 
with brightness smaller than that of the star by (R p /a) 2 . 

The effect of the crescent that we calculated, which result in shoulders around transits, 
is therefore not offsetting the light refracted away during transit. The latter is offset mostly 
by deeply penetrating rays, that would allow that planet to shine anywhere in the orbit. The 
shoulders provide only a small correction to the above. Together, averaged over all possible 
planet positions, the shoulders and the deeply penetrating rays which show up away from 
transit, cancel the decrease of light during transit exactly. However, in any realistic situation 



This average is not over an orbit, but over all orbits with all inclinations 
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deeply penetrating rays will be absorbed, and therefore less light, on average, will arrive to 
the observer. 



5. Extinction by Rayleigh scattering 

In order to determine the relevance of refraction one must estimate the starlight extinc- 
tion as it passes through the planet's atmosphere. Rayleigh scattering is a minimal cause of 
extinction. We ignore molecular absorption and clouds. We also ignore limb darkening, and 
denote the uniform, unextincted, stellar surface brightness by Jo- The observed brightness 
after a given path through the planet's atmosphere is / = Joexp (— r), where r is the optical 
depth: 

T = J aNdl. (31) 

Here, a = f^r^^r- is the total cross section for Rayleigh scattering per molecule of the 
gas, n is the index of refraction, iV is the number of molecules per unit volume, A is the wave 
length. The integration takes part over the path of light. Here we assume that the planet 
has a Hydrogen atmosphere, in this case we have 

d = ^1 = 3.046^ (32) 
N mol 

For an atmosphere in hydrostatic equilibrium the gas density decays exponentially over the 
scale height H <C R p , so that the principal contribution to the optical depth comes from the 
integration over a few scale heights near the impact parameter b = R p . Therefore, for a light 
ray with negligible deflection, 

r M ( ypT¥-b* \ 

r = J aN max exp I — I dl, (33) 

where / is the coordinate along the light trajectory measured from the point of closest 
approach to the planet center. By expanding the root around small l/b, and assuming 
b ^> H, we obtain an analytic expression for the optical depth: 

cy . 32n 3 d TT . ^ . lS 

T = ^ Ae = -Jv HAe (34) 

Now, using the expression for the cross section o and assuming that the planet's tem- 
perature is determined from thermal equilibrium with albedo A, 

T = g^MWV /2 (1 - XrAe . (35 ) 



3V2X 4 iiGM p V a 
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We now derive asymptotic results far away from transit, X ^> 1, as well as just before 
transits R p /R* < J-l< 1. 



5.1. Away from occultation regime 

When I>1, the refraction angle of the close and distant edges of the star do not vary 
considerably, AO = Therefore, from equation (1351) . the optical depth is uniform. The 
crescent specific intensity is constant, attenuated simply by e~ T , where 

0.067(VV^V /4 ( Rp 



A / \ Lq J \ i?j 

3/2 

X 



Here we have used red light with A re d = 650 nm as our fiducial wavelength. Clearly, the 
optical depth in the blue part of the spectrum is significantly larger. 



The fractional flux from the crescent before extinction, is given by equation (1231) . it 
decreases with the planet's semimajor axis a as a _1//2 . On the other hand r decreases as 
fl -3/2 t ransm j^ ec l light is therefore maximal when 



We therefore define 



a 1/3 = 0.34 AU 



1/3. (37) 



8/3 / r \ 1/6 / p \ 4/3 



^red \ * \ / n P 



Ay \ Lq J \ R) 



to be the semimajor axis where optical depth is 1/3 for deflections at angles i?*/a. For our 
fiducial parameters, ai/3 corresponding to an orbital period of 74 days. 

For a > the refraction effect will be practically unextincted, with magnitude given 
by equation (123]) . for any X < X max given by 

X max = (a/a 1/3 ) 2/3 . (39) 
For X > X max the crescent flux will decay exponentially with X due to Rayleigh scattering. 
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5.2. Near occultation regime 

When the planet center as projected on the sky is less than a stellar radius away from 
the star's edge (point A in figured]), the refraction angles of the close and distant edges of 
the star (points A and C on figure [1]) differ considerably. The light from point A is deflected 
much less than R^/a and therefore will be significantly less extincted than we estimate in 
§5.11 Hence, in order to estimate the attenuation of the star light that is refracted through 
the planet atmosphere, one needs to integrate along the crescent. In that case the fractional 
size of the effect, taking extinction into account, is given by 



J 



S'J^e-^dr _ Rp /■« 



S f^dr Rl 



out 



T(r W, (40) 

-Win 



'-Rin 

where R out and R{ n the maximal and minimal radius of the crescent as measured from the 
center of the planet, respectively. The denominator integral is simply the crescent radial 
width in the g < X - K 1 regime. Changing variables from r to r we obtain 
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Here T in and r ou t refer to the optical depth associated with the refracted light that sets the 
inner and outer edges of the crescent. Explicitly 

r m = ^HA& A , T out = ^HAQ c > (42) 

From equation ( |27|) we see that r in is given by equation ( 1361) with X = 2 while r out is obtained 
if X is replaced by X — 1 + R p /R*. We examine the attenuation scaling in a few different 
regimes. 

• Both Tin, Tout 1- In this regime, the crescent luminosity will undergo considerable 
extinction: 

S' HR„ e~ T ° ut . . 

— = £ (43) 

S ext R% T out 

This applies for a {2R p /?)R jf ) 2 ^ai/- i or, for our fiducial parameters a 0.06 AU, 
i.e., a period shorter than 5 days. Most current transits, therefore, fall in this category, 
making refraction effects difficult to see. The additional light in this case grows expo- 
nentially in time as the planet approaches transit. Most light arrives when the planet 
is about R p away from transit. Therefore, the duration of the increased light would be 
shorter than the duration of the transit by R P /2R*, i.e. about a factor of 20. For hot 
Jupiters, transits, which last about 3 hours, should be accompanied by slightly bright 
edges for about 10 minutes. 
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• The crescent is only partially extincted. r m ^> 1 but r out <C 1. Here the crescent 
luminosity changes considerably along its radial width, so the crescent luminosity is 
dominated by regions where r < 1. We obtain: 

^ ~ -ff-Rp j / 1 | (44) 

S ext R+ \T ou t J 

relevant for 3~ 2 ^ 3 a 1 / 3 ^> a^> (2R P / '3-R*) 2, /3 ''ays or for our fiducial parameters 0.16 AU > 
a > 0.06 AU or periods between 5 and 25 days. The fractional flux undergoes two com- 
peting effects. On one hand the area of the crescent decreases as a _1//2 , on the other 
hand r out oc a _3//2 decreases with a, allowing a larger fraction of the light to pass 
through. The maximum is obtained where r out = e~ 3 = 0.05. However, at such large 
semimajor axis, given our fiducial parameter, T- m = {R*/R p )t ou ± = 0.5 < 1. There- 
fore, the maximum fractional luminosity of the crescent is obtained at the transition 
between this and the next regime, when r in = 1, with orbital periods of 25 days. There, 

j = ^ HRJR P ) = 3.6 x 10- 5 . (45) 

This is the largest effect refraction can produce. The numerical value can be slightly 
higher if the planet's radius is larger than that of Jupiter and if its mass is somewhat 
lower, allowing for a larger scale height at a given temperature. 

• Both ri n , r out ^ 1- Here the crescent luminosity is hardly attenuated, so that its 
fractional luminosity is given by ( 13"U|) . It decreases with semimajor axis showing again, 
that the maximum is obtained around r in = 1. 



6. Summary 

We have derived the possible paths of light passing through a planet with exponential 
deviation of the index of refraction from unity. Though we explore exotic trapped rays and 
rays following circular orbits, we focus on small deflections, as given by our equation f[T2"j) . 

We explore the size and shape of the image of the star as refracted through the planet's 
atmosphere. The size of the image is given by equation ( 123]) far away from transit while 
the largest and least extinct effect should be seen just before transit, and its magnitude is 
given by equation (I3"U|) . For planets with short orbital periods, refraction involves relatively 
large deflections, and will therefore be reduced by extinction. On the other hand, planets 
with large orbital period have cold surfaces resulting in a small effect. The most noticeable 
effect would be for planets with orbital period of about 25 days. For extrasolar planets with 
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R p = 1.5Rj, like HD 209458 where R p = lARj and M = 0.7Mj (Brown et all 120011 ) . the 
fractional flux increase would be as large as 10~ 4 , if it was located on a 25 days orbit. This 
is readily observable with current instruments. 



Hui fc Seagerl (120021 ) have also considered refraction by planetary atmospheres. How- 
ever, they focused on caustics, which are important for point sources. These are relevant for 
refraction of background stars by atmospheres of planets in our solar system. We instead, 
were interested in the case where the source is very large, smearing out the effects of caustics. 
We have also ignored the effect of oblateness that they considered. 

For HD 209458 on its observed 3.5 days orbit, the refraction effect would be extincted 
with optical depth of about r out = 1.7. The refractive effect, even given this extinction would 
be about 10 -5 . It would last for about 10 minutes just before and just after transit. This is 
still detectable in existing or future data. If the enhancement of light is detected, it could 
provide a direct measurement of the scale height of extrasolar planet. 

The timescale for the brightest part of the refractive effect is somewhat shorter than 
the duration of the transit, especially if extinction is involved for planets with short orbital 
periods. 

We also show, that during the eclipse, refraction of stellar light away from the observer 
may cause a signal similar to that of a regular geometric eclipse. We define the radius R p as 
the effective radius of the planet as observed by a refractive eclipse. Indeed, for any planet 
on orbital periods longer than 25 days, R p , the radius where refraction effect control the 
transit occur higher up in the atmosphere that R ray , where Rayleigh scattering is important. 
Using equation (1381) we find that these two radii are related by 

R p = Rray + ^H\n( — ) (46) 
2 V a V3 / 

Of course, other kinds of extinction may be in place, most notably clouds, which we have 
ignored. Those may increase the effective extinction radius above R ray and perhaps beyond 
our effective radius R p . 

Finally, we mention that refraction effects would have only a slight modification on the 
Rossiter - McLaughlin effect, but perhaps a large effect on the absorbed spectrum during 
transits. If R p — R Tety ^> H such that refraction rather than extinction determine the eclipse, 
less absorption would be observed in the spectrum. 
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